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PACS 46 . 55 . +d - Tribology and mechanical contacts 
PACS 81.40.Pq - Friction, lubrication, and wear 

PACS 85.85.+j - Micro- and nano-electromechanical systems (MEMS/NEMS) and devices 

Abstract. - We report on normal stress field measurements at the multicontact interface between 
a rough elastomeric film and a smooth glass sphere under normal load, using an original MEMS- 
based stress sensing device. These measurements are compared to Finite Elements Method cal- 
culations with boundary conditions obeying locally Amontons' rigid-plastic-like friction law with 
a uniform friction coefficient. In dry contact conditions, significant deviations are observed which 
decrease with increasing load. In lubricated conditions, the measured profile recovers almost per- 
fectly the predicted profile. These results are interpreted as a consequence of the finite compliance 
of the multicontact interface, a mechanism which is not taken into account in Amontons' law. 



Knowledge of the stress field at the contact region be- 
tween two solids is of considerable interest to numerous 
domains such as mechanical engineering [1-6], solid fric- 
tion [7-9] or seismology [10]. From a continuum mechanics 
point of view, the theoretical or numerical determination 
df this field requires a set of constitutive equations char- 
acterizing the mechanical response of the interface. These 
are usually inferred from macroscopic measurements: for 
instance, frictional contacts are often found to obey the 
empirical Amontons' friction law which states that irre- 
versible sliding occurs when the ratio of tangential to 
normal forces reaches a static friction coefficient jimacro 
without any prior deformation of the interface [11, 12]. 
The classical approach for calculating frictional contact 
stress field therefore consists in considering a smooth in- 
terface exhibiting an analog rigid-plastic response: fimacro 
now defines the threshold ratio between shear and normal 
stress components for local slip to occur [5,7-9]. 

Considering both a smooth interface and Amontons' law 
may seem paradoxical since the latter is expected to be 
valid only for rough solids. Because of adhesion forces, 
contact between molecularly smooth surfaces generally do 
not obey Amontons' friction law [13, 14]. Furthermore, 
random roughness provides a microscopic basis to Amon- 



tons' law [15-18]: it allows for the description of the inter- 
face as a collection of isolated load-bearing points forming 
a multicontact interface (MCI). For standard roughness 
characteristics, the real area of contact - and thus the tan- 
gential force required to trigger sliding - grows linearly 
with the applied load. 

When associated with Amontons' law, the smooth in- 
terface hypothesis must therefore be understood as the 
limiting case of a macroscopic surface bearing vanishingly 
small roughness. One may wonder to what extent this ap- 
proximation is valid in real MCI's. Experimentally, this 
question has been addressed by focusing on the relation- 
ship between the macroscopic normal load and the inden- 
tation depth [17,19-21]. Here we go further and compare 
the stress field measured within an extended MCI with Fi- 
nite Elements Method (FEM) calculations assuming both 
smooth surfaces and the local (rigid-plastic) Amontons' 
law. 

In order to perform such local measurements, we take 
advantage of the recent development of Micro Electro Me- 
chanical Systems (MEMS) [22,23]. Local normal stress 
measurements are obtainend using a MEMS force sen- 
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Fig. 1: Schematics of the stress sensor. A MEMS force sen- 
sor enables to measure the applied normal stress at the base 
of a rough, nominally flat PDMS film (thickness h — 2 mm, 
lateral dimensions 50 x 50 mm) pressed against a spherical 
glass substrate. The macroscopic normal and tangential loads 
are measured through the extension of two orthogonal loading 
cantilevers (normal stiffness 641 ±5 iV.ra -1 , tangential stiffness 
5 1100 ±700 iV.ra -1 ) by capacitive position sensors (respectively 
MCC30 and MCC5, Fogale nanotech). 



sor 1 embedded at the rigid base of a rough, nominally 
flat elastomeric film pressed against a rigid sphere under 
normal load (Fig.l). The MEMS sensitive part is a rigid 
cylinder (diameter 550 //m, length 475 jam) attached to a 
suspended circular Silicon membrane (radius 1 mm, thick- 
ness 100 jam) whose deformations are measured with four 
couples of piezo-resistive gauges (see inset of fig.l). The 
applied normal stress is therefore averaged over a surface 
of a few millimeters square. 

The elastomeric material is a cross-linked polydimethyl- 
siloxane (PDMS) (Sylgard 184, Dow Corning) of Young's 
modulus 2.2 ± 0.1 MPa, and of Poisson ratio 0.5 [25]. No 
measurable stress relaxation being observed after a sudden 
loading, the PDMS can be considered as purely elastic. 
The film is obtained by pouring the cross-linker/PDMS 
melt on the MEMS into a parallelepipedic mold covered 
with a PolyMethylMethAcrylate plate roughened by abra- 
sion with an aqueous solution of Silicon Carbide pow- 
der (mean diameter of the grains 37/xra). After curing 
and demoulding, the resulting rms surface roughness is 
measured with an interferential profilometer (M3D, Fo- 
gale Nanotech) to be p = 1.82 ± 0.10 jim. This roughness 
is sufficient to avoid any measurable pull-off force against 
smooth glass substrates, as discussed in [26]. 

The stress sensing device is calibrated by indenting the 
film surface with a rigid rod of diameter 500 /im, under 
a normal load P. With this type of indentor, the sensor 
output is found to be linear with the applied load. By 
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Fig. 2: a) Radial profiles associated to the indentation of the 
elastic film by a rigid thin rod, under normal load P = 1 N. 
The measured output voltage U zz (•) is compared to FEM 
results (solid line) for the normal stress a zz at the base of the 
elastic film. The error bars represent the electronic noise. The 
black rectangular patch represents the rod diameter (500 fim). 
b) Normalized apparatus function of the MEMS for the normal 
stress, G zz (•), as determined by numerical integration of eq.l, 
and its Gaussian fit in dashed line. 



successively varying the position of the rod along the x 
direction, and assuming homogeneity of the surface prop- 
erties of the film, one can construct point by point the 
normal output voltage radial profile U zz (fig. 2a)). The 
latter is then compared to the results of Finite Elements 
Method (FEM, Software Castem 2007) calculations for the 
stress a zz at the base of a smooth axi-symmetrical elas- 
tic film (of same elastic moduli and thickness) perfectly 
adhering to its rigid base and submitted to a prescribed 
normal displacement over a central circular area 2 of diame- 
ter 500 fim. As expected for contact regions of dimensions 
smaller than the film thickness, the stress calculated at the 
base of the film is found to be insensitive to the frictional 
boundary conditions. 

Ignoring the stress field modifications induced by the 
MEMS 3D structure, one can relate the measured output 
voltage U zz to the stress field a zz by writing down that 

U zz (x, y) = A ZZ G ZZ <g> a zz (x, y) (1) 

where A zz is a conversion constant (expressed in mV/ Pa), 



1 This MEMS also measures the tangential stress, but with an 
accuracy which is not sufficient for the present study [24]. 



2 Such results could have been obtained semi-analytically for fric- 
tionless conditions by using the model developed in [27] but FEM 
calculations have been preferred here because they allow for variable 
boundary conditions. 
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G zz is a normalized apparatus function and is a convo- 
lution product. In Fourier space, eq.l becomes 



A zz G zz (x,y) = T 



F{Uzz}(fxJy) 
F{°zz} (fx,fy) 



(x,v) ( 2 ) 



where T is the bidimensional spatial Fourier Transform, 
T~ x its inverse, and where f x , f y are respectively the spa- 
tial frequencies in the x, y directions. The U zz (x,y) and 
a zz (x, y) fields are built from the profiles along the x axis, 
assuming axisymmetry, and then transformed using the 
Fast Fourier Transform (FFT) algorithm. The fast decay- 
ing of T {<J zz } with increasing spatial frequency introduces 
divergences of the ratio in eq.2. To avoid it, a white noise 
of amplitude ten times weaker than the weakest relevant 
spectral component is added to both terms of the ratio 
before applying the FFT. The result is found to be insen- 
sitive to the particular amplitude of this white noise. A zz 
is determined so that the integral of G zz is 1. 

Figure 2b) shows the resulting apparatus function. It 
has a bell shape with typical width of the order of 600 fim 
comparable to the MEMS lateral dimension. For the sub- 
sequent convolutions an approximated apparatus func- 
tion, taken as a gaussian of standard deviation 561 fim 
(fig. 2b)), has been used. This approximated apparatus 
function gives back the measured profile when convoluted 
by the calculated one, proving that neither the deconvolu- 
tion method nor the approximation introduces significant 
loss of information. In particular, the point-like indenta- 
tion involves a large enough spatial spectrum to allow for 
a faithfull determination of G zz within all of its relevant 
spatial components. 

Sphere-on-plane MCIs are formed against an optical 
plano-convex spherical glass lens (radius of curvature 
128.8 mm). Both the glass and PDMS surfaces are pas- 
sivated using a vapor-phase silanization procedure which 
lowers and homogenizes the surface energy. Each contact 
is obtained using the following loading sequence. The glass 
lens is pressed against the film up to the prescribed load P 
within 2 % relative error. Due to this loading, the extrem- 
ity of the normal cantilever is tangentially displaced and 
a significant tangential load Q is induced. Consequently, 
from this position, the contact is renewed by manual sep- 
aration which results in a much smaller but finite Q. The 
glass lens is eventually translated a few micrometers tan- 
gentially down to Q = 0. Both the surface treatment 
and the loading sequence are found to yield an excellent 
reproducibility of the measurements. As for the rod in- 
dentation, the radial profiles are derived from a series of 
33 contacts whose centers lie every 0.5 mm along the x di- 
rection. These profiles divided by A zz have the dimension 
of a stress and are labelled S zz . 

For a quantitative comparison, FEM calculations are 
carried out for a frictional sphere-on-plane contact with 
the same geometry. Both contacting surfaces are taken 
as smooth and the interface is assumed to obey locally 
Amontons' friction law with a friction coefficient \i. Both 



S , G <g)<7 (kPa) 

ZZ ZZ ZZ a \ 





x >jp 








-2 2 

x {mm) 



Fig. 3: Measured normal stress profiles (S zz ) under normal 
loading by a rigid sphere (P = 0.34 N (•), 1.37 N (■), 2.75 N 
(A)). Comparison is made with G zz ® a zz for \i — (solid 
lines, indentation depths 18, 33 and 45 fim respectively) and 
H = oo (dashed lines, indentation depths 16, 28 and 37 \ivn 
respectively), a) Linear scale, b) Semi- logarithmic scale. The 
black rectangular patches on a) represent the contact diame- 
ters (2.00, 2.90 and 3.45 mm for P = 0.34, 1.37 and 2.75 N 
respectively) obtained from the FEM calculations for /x = 0. 



solids are discretized with a uniform mesh size of 50 jim 
and the normal displacement of the rigid elastic sphere 
is prescribed. The contact conditions are satisfied using 
a double Lagrange multiplier implying that both surfaces 
are slave and master. The normal load is reached step- 
by-step and at each step an iterative Newton-Raphson 
method is used to satisfy both the unilateral contact and 
the friction law. 

Figure 3 compares the S zz and G zz ®a zz profiles for two 
limiting boundary conditions, fi = and fi = oo, and for 
three values of P. Within the error bar the measured pro- 
files are bracketed by the two limiting numerical profiles 
over the whole spatial range and over 3 orders of magni- 
tude, as clearly displayed on fig. 3b). In the contact outer 
region, S zz is systematically very close to the frictionless 
profile whereas at the center it increasingly departs from 
it with the load, as discussed further. Similar measure- 
ments are performed, for two limiting loads (P = 0.69 N 
and P = 2.75 A 7 "), under lubricated conditions. A glyc- 
erol droplet is inserted at the interface prior to loading. 
The profiles display both a larger spatial extent and a 
lower maximum amplitude than the ones obtained with 
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Fig. 4: Measured normal stress profile (S zz ) under normal 
loading (P = 2.75 N) by a rigid sphere for both dry (•) and 
glycerol- lubricated (■) contacts. Shown in solid and dashed 
lines are the G zz a zz profiles for 4 values of \i (0, 0.3, 0.4, 
oo). The black rectangular patch represents the contact diam- 
eter (3.45 mm) obtained from the FEM calculation for \i = 
(indentation depth 45 /am). 

dry conditions. They are correctly captured by the FEM 
calculation using a null friction coefficient (as shown for 
P = 2.75 7V on fig.4). 

Except for the smallest load, any pressure profile under 
dry conditions can be correctly approached using an effec- 
tive friction coefficient p e which can be determined by trial 
and error with a resolution of order 0.1 (as illustrated in 
fig.4). The load dependence of p e is exhibited by plotting 
the maximum pressure S max normalized by the maximum 
pressure p max calculated for a frictionless contact, for dif- 
ferent loads. On the same graph, we plot the equivalent 
expression p max /p max obtained by FEM calculation for 
different friction coefficients p in the range [0, oo]. This 
representation allows one to directly read the value of the 
effective friction coefficient p e for any load. p e is found to 
increase with P but remains well below the macroscopic 
friction coefficient p ma cro ~ 1.6 over the range of loads 
explored 3 . 

The significant discrepancy between the effective and 
the macroscopic friction coefficient suggests that the rigid- 
plastic response of the interface underlying Amontons' 
law needs to be refined. One may in particular question 
the rigid hypothesis since MCI's are known to possess fi- 
nite compliances both in normal and in-plane directions. 
The effect of the normal compressibility of a rough inter- 
face [16,17,21,29] on a sphere-on-plane contact has first 
been described by Greenwood and Tripp [30]. They pre- 
dict an increase of the apparent contact radius with re- 
spect to a smooth interface, as well as a decrease in the 
maximum normal stress. These deviations are expected 
to vanish when the ratio of the rms surface roughness p 

3 the macroscopic friction coefficient has been measured for a driv- 
ing velocity v = 100/ira/s of the rigid base of the sensor: it shows a 
(small) decrease from 1.8 to 1.5 as the load is increased from 0.34 N 
to 2.75 N, a behavior usually attributed to the finite adhesion energy 
of the interface [28] 




Fig. 5: Load dependence of Smax, the maximum normal stress 
measured for dry (•) and lubricated (■, guideline in dashed 
line) conditions, and p max , the maximum calculated normal 
stress for various values of p. Both quantities have been nor- 
malized by p maX i the value of p max for p = 0. The vertical 
arrows represent the apparent drop of the friction coefficient 
from \imacro ~ 1.6 down to p e for dry conditions. 

to the indentation depth becomes small. The effect of the 
MCI tangential compliance has been recently probed ex- 
perimentally [31,32] : a global reversible deformation of 
the interface between two contacting solids is measured 
before irreversible slippage occurs ; regardless of the nor- 
mal load, the maximum shear deformation of the interface 
before slippage is of the order of its rms roughness p. 

In the present set of experiments, these two effects can 
be uncoupled since the tangential stress at the interface 
vanishes under lubricated conditions. The fact that the 
lubricated profiles (■ on fig. 5) are compatible with those 
calculated for a frictionless contact suggests that the MCI 
normal compressibility effect is negligible in most of the 
load range explored (typically for P > IN). Still, it is 
probably responsible for the fact that Smax/ Pmax ^ a ^ s be- 
low 1 for the smallest loads when the indentation depth 
becomes comparable to the thickness p of the rough layer. 

The tangential compliance of the interface is thus ex- 
pected to be primarily responsible for the observed drop 
from Prnacro to /i e for dry contacts. This effect can be 
qualitatively understood by first considering the case of 
an interface with infinite friction. The normal loading of 
a sphere-on-plane contact yields a divergence of the shear 
stress at the edge of the contact region [5,6]. In a system 
with finite friction, this stress is relaxed by the develop- 
ment of a slip annulus at the periphery of the contact re- 
gion which coexists with a central (circular) stick region, 
as discussed by Spence [5] and observed in our FEM calcu- 
lations. The radius of the stick region is a growing function 
of the friction coefficient. As compared to a rigid interface, 
the existence of a finite tangential compliance allows for 
a subsequent strain relaxation at the elastomeric surface 
which is qualitatively equivalent to reducing the value of 
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the apparent friction coefficient. This effect should van- 
ish as the ratio of the tangential displacement induced by 
the sphere-on-plane loading becomes larger than the max- 
imum elastic displacement p allowed by the MCI defor- 
mation. The effective friction coefficient is thus expected 
to asymptotically reach pmacro as the load P is increased 
in qualitative agreement with our observations (fig. 5). It 
is noticeable however that the reduction of the apparent 
friction coefficient remains significant even for loads asso- 
ciated with a normal displacement of the rigid substrate 
more than 20 times larger than the interface thickness p. 

The results reported in this Letter provide evidence for 
the inadequacy of local Amontons' friction law to correctly 
capture sphere-on-plane MCI pressure profiles. They pro- 
vide the first quantitative experimental measurements of 
such stress fields over a large range of loads for both lu- 
bricated and frictional contacts. It appears that each pro- 
file can be correctly described within the scope of Amon- 
tons' rigid-plastic law, but with a load-dependent effective 
friction coefficient. This coefficient grows with the total 
applied load but remains much smaller than the macro- 
scopic friction coefficient even when the standard smooth 
hypothesis criterion (rms surface roughness much smaller 
than the indentation depth) is reached. By comparing lu- 
bricated and dry frictional contacts, we can separate the 
effects of normal and tangential compliance of the MCI 
and conclude that the latter is primarily responsible for 
the observed deviation to the local Amontons' law. 

More generally, this study provides experimental sup- 
port to previous works aiming at taking into account 
the rheological behaviour of the micro-contacts popula- 
tion in the mechanical studies of multicontact interfaces 
under tangential load [33-35]. From the present results 
we suggest that it should also be done for contacts un- 
der purely normal load. A first step could be the use of 
an elasto-plastic-like friction law instead of the classical 
rigid-plastic-like Amontons's one. These results were ob- 
tained with an original MEMS-based stress sensor, which 
has proven to be well-suited for the study of stress fields 
in centimeter-sized contacts and could be used to test di- 
rectly any other mechanical model of the frictional inter- 
face. Many other aspects of contact mechanics are likely 
to be probed with the same method, such as the dynami- 
cal frictional regimes or the history-dependence of a con- 
tact submitted to an oscillatory tangential load below the 
sliding threshold. Other domains such as rheology or ad- 
hesion, where accurate spatially resolved stress measure- 
ments at interfaces are needed, could also benefit from the 
approach described here. 
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